Creating and observing TV-partite entanglement with atoms 
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The Mermin inequality provides a criterion for experimentally ruling out local-realistic descrip- 
tions of multiparticle systems. A violation of this inequality means that the particles must be 
entangled, but does not, in general, indicate whether A-partite entanglement is present. For this, 
a stricter bound is required. Here we discuss this bound and use it to propose a scheme for demon- 
strating A-partite entanglement with Rydberg atoms in microwave cavities. 
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Entanglement is a central feature of quantum mechan- 
ics and the key resource in a wide range of quantum infor- 
mation processing tasks. Being able to detect and charac- 
terize the entanglement present in a system is therefore 
an important challenge for physicists. The first math- 
ematically sharp method of detecting entanglement for 
pairs of particles was proposed by John Bell in 1964 []]]. 
His now-famous inequality provided an unambiguous way 
of distinguishing quantum-mechanical predictions from 
those of local realistic models. Multiparticle generaliza- 
tions of the Bell inequality were subsequently provided by 
Mermin Q and others 0, 0, 111 ■ These can be used to rule 
out local-realistic models for A^-particle systems and are 
also interesting because there is a close relationship be- 
tween their violation and the security of A-partner quan- 
tum communications [f|. 

In general, a violation of these inequalities means that 
the particles must be entangled. In fact, by placing a 
stricter bound on the inequality 0,1! it is even pos- 
sible to determine what class of entanglement is present 
from 2-entangled to A-entangled states. This method 
has been used to experimentally confirm three-body en- 
tanglement for photons [Icj . In this Letter we discuss 
the conditions required to demonstrate iV-partite entan- 
glement and show how it could be generated and detected 
in atomic systems. We begin by reviewing the Mermin 
inequality and discussing how A^-partite entangled states 
could be generated with Rydberg atoms in microwave 
cavities. We then demonstrate how all the terms in the 
Mermin inequality could be measured and consider some 
of the practical issues surrounding such an implementa- 
tion. 

It is helpful to start with a brief overview of the Mer- 
min inequality 0. For this, we consider an A^-particle 
GHZ state of the form [Hj], 



state, |$)jv, is an eigenstate of the operator, 
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where g or e in the jth position denotes that the jth par- 
ticle is in the ground or excited state. In an experiment, 
these Af particles would be spatially separated and mea- 
surements made on each of them. We note that the GHZ 
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3 = 1 



(2) 



with eigenvalue 2 N ~ 1 , where o l x and a y are respectively 
the x and y Pauli spin matrices acting on the jth. particle. 
Expanding, F = {^[An]^), for some general A-particle 
state l^), one finds, 



F= 1 - 
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where each line of ([3]) contains all distinct permutations 
of the operators in the term shown on that line. Only 
terms with an even number of a y operators are non-zero. 
For the GHZ state, we have F = N {$\A N \$) N = 2 N ~ 1 . 

We would now like to find the corresponding value of 
F for a local hidden-variable state. Following Mermin, 
we consider the case where the measured distribution 
functions, P^...^ (mi • • • ttin), where £ {x,y} and 
rrij e {<7, e}, that describes the 2 N ~ l measurements that 
must be performed to yield the correlations in ([3]), can 
be written in the conditionally independent form, 



P.. 



(mi • • • m N ) = / dXp(X) [p* x (mi , A) • 
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where (mj, A) is the probability distribution of results 
for a measurement in the /ii basis on particle i. This 
general hidden-variable form attributes the correlations 
to some unspecified set of parameters, A, common to all 
N particles, with distribution p(X). It accounts for cor- 
relations in terms of information jointly available to the 
particles when they left their common source. An anal- 
ysis similar to that of Mermin's shows that such a local 
hidden- variable state has the bounds 



F < 2 N/2 , N even, 

F < 2< Ar - 1 )/ 2 , iVodd. 



(5) 
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We see that the GHZ state, \<&) N , for which F = 2 N ~ 1 , 
violates ([5]) by an exponential amount. 

A violation of inequality ([5]) demonstrates that the par- 
ticles are entangled, but does not guarantee TV-partite 
entanglement For this, a tighter bound is needed, 

which can be found as follows. We assume that there 
is, at most, (N — l)-partite entanglement in the system. 
This means that the state of one of the lattice sites (say 
the first one) factorizes and the expectation value of the 
corresponding operator factorizes in each term in (|3j). 
Then using the fact that all expectation values disap- 
pear if they contain an odd number of a y operators, it is 
straightforward to show that we are left with: 



F = 



< 2 



N-2 



(6) 



where the last step follows because (a x ) < 1 and Fn-i < 
2 N ~ 2 . This agrees with the condition found for three 
particles in |8j. Any state that violates ([6]) must have at 
least A-partite entanglement. 

The aim of an experiment would be to measure all the 
expectation values in ([3|) and see firstly whether their 
sum violates the bound given by (O, in which case en- 
tanglement is present, and secondly whether it violates 
the bound given by (0 in which case A-partite entan- 
glement is present. We now discuss how such a scheme 
could be implemented. 

The first step is to create a GHZ of the form of (JTJ) . 
To do this we use a modified version of the scheme found 
in [ijj]- In short we start with A atoms of 85 Rb, each 
initially in the Rydberg state 63P. Using resonant mi- 
crowave fields transitions can be driven to the levels 6ID 
and 62P. We shall refer to 63P as |e), the relative excited 
state, and 61D as \g), the relative ground state. These 
form the basis used to construct the GHZ state. The 
state 62P is referred to as \i), an auxiliary state. The 
state |e) is produced using a three step laser excitation. 
Zheng and Guo demonstrate in their paper [l3| that us- 
ing three states like this, with a cavity detuned from the 
|e) <-* \g) transition it is possible to create an EPR state. 
This was later demonstrated experimentally by Osnaghi 
et al. [3]. Recently this scheme was extended in [l2[ to 
show that using the same states and A — I cavities in a 
line it is possible to produce a GHZ state of N atoms [23] • 
An alternative has been demonstrated by Rauschenbeu- 
tel et al. 



15] 



The scheme in [12j starts with all atoms initiall y in 
the state |e). These atoms are produced on demand [16j |. 
Next all atoms are rotated to the state |+) = (\g) + 
\e))/^/2 using a microwave field resonant with the |e) <-> 
\g) transition. This corresponds to a rotation about the 
x-axis on the Bloch sphere. The Hamiltonian used to 
implement this is 



Hi = nn (<t+ + &-) 



(7) 



where ft is the coupling strength of the atom with the 
field and a + and a~ are the atomic raising and lowering 



operators. The amplitudes of the atom given as a(t) \g) - 
b(t) \e) evolve according to the equations 



a(i) = clq cos(fii) — ibo sm(Qt) 
b(t) = bo cos(fit) — iao sin(Slf) . 



(8) 



which is equivalent to a rotation about the x-axis. 

The first atom, which will interact with each other 
atom in turn, then passes through another microwave 
field (labelled as b in Fig. [T|) which is resonant with the 
transition |e) <-*■ \i) leaving the first atom in the state 
(\g) + \i))/\/2. The first atom now interacts with the sec- 
ond in a high-Q microwave cavity. This cavity is detuned 
from the |e) <-> \g) resonance. The four possible states 
that may interact are shown in the left column of table 
HJ The first two states in the left column of the table are 
modelled using a two atom Tavis-Cummings model with 
large detuning. With a zero photon field [lTf this can be 
solved to give the interaction 



|S, e) i-> e ' nt [cos(7i) \g, e) - i sm(<yt) \e, g)} 

\g,g) ^ \g,g) 



(9) 



where 7 = g 2 /A, g is the atom-field coupling constant 
and A is the detuning. As the field is detuned from the 
transition, excited atoms can only virtually excite the 
field, effectively coupling the atoms together and allow- 
ing an excitation to be passed between atoms. The zero 
photon stipulation is satisfied in the laboratory by cryo- 
genically cooling the cavity, and possibly by preceding 
the experiment with a chain of atoms in the state \g) 
to unload the field. The second two states in the left 
column of table U each have one atom in the state \i), 
for which all transitions are so far from the cavity reso- 
nance that we may assume that it does not contribute to 
the dynamics of the system. This leaves one atom and 
the detuned field, which is modelled using the Jaynes- 
Cummings model. In the case of zero photons and large 
detuning the evolution of these two states is given by 



\i,e) 

\i,g) \i,g) , 



(10) 



For creating a GHZ state we choose the interaction time 
and detuning such that t — tt/Aj. For each interac- 
tion this provides the right hand side of the truth table 
HI After the interaction with the first atom each atom 
emerging from each cavity is rotated so that \g) — > |+) 
and |e) -> ]-) = (\g) - |e))/\/2. After the first atom 
emerges from the final cavity it enters a microwave field 
that drives the transition |z) <-» |e) so that the atom is 
left with no amplitude in \i). The atoms are now in a 
GHZ state ready for use. Figure Q] demonstrates the con- 
struction of a three atom GHZ state. It is also possible to 
generate larger entangled structures with more elaborate 
cavity arrays [l~8l |. 
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Input Output 

\g,g) \g,g) 

\g,e) \g,e) 

\i,g) \i,g) 

\i,e) -\i,e) 

TABLE I: The truth table for the collisional phase gate in- 
troduced in 
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FIG. 1: This example produces a three party GHZ state of 
atoms. Atoms approach as indicated by arrows to coincide 
in each cavity (large circles). Initially all atoms are in state 
je). The rotation zones labelled a take \e) to |+) and c take 
|e) to |— ) and \g) to |+). Those zones labelled b act on the 
first atom (traveling horizontally) to switch \e) components 
of the state to \i), the auxiliary state, and vice versa. The 
enclosed region is where the GHZ state exists. The zones la- 
belled d are arbitrary rotations used to study the state before 
measurement. 

Now we would like to make measurements on this sys- 
tem that correspond to the terms in ([3]). This involves 
making measurements on each qubit in the basis of the 
eigenstates of a x and a y , i.e. |g)±|e) and \g)H |e) respec- 
tively. The rotations take place in the regions labelled d 
in Figure [TJ The rotation in the basis of a x has already 
been discussed, and only one other rotation on the Bloch 
sphere is needed to realise a rotation in the basis of a y . 

The second operation we use is a rotation about the 
z-axis of the Bloch sphere, R 2 . This can be implemented 
simply by applying an electric field to the atom. This 
effectively increases the transition energy between the 
levels of the atom and results in a different phase evo- 
lution for the two states. This technique has been exper- 
imentally demonstrated with sodium Rydberg atoms by 
Ryabtsev et al. [HJ . They applied a resonant microwave 
pulse to perform an x-rotation between two sodium Ry- 
dberg atoms followed by a Stark shift z-rotation and an- 
other x-rotation. This combination of interactions al- 
lowed them to perform Ramsey interferometry of the R 2 
operation applied by the Stark shift. Their results con- 
firmed that they has successfully implemented the phase 
operation R 2 and that it is coherent [l9[ ■ 



For Rydberg states of alkali metals the Stark shift 
varies approximately quadratically with the electric field 
strength [H, 

Sg ie (X -Ctg^E 2 , (II) 

where a 9te is the polarizability of the ground or excited 
state of the atom and E is the amplitude of the electric 
field. The polarizability of the two states will be different, 
so the relative phase shift will be given by 

9 = e(a g - a e )E 2 , (12) 

where e is a factor calculated by integrating over the pulse 
shape of the electric field. Applying the electric field for 
some time t, the state of the atom evolves as 

a\e)+b\g)^e-* et a\e)+b\g). (13) 

This corresponds to a rotation about the z-axis, where 
<dt is the angle of rotation. 

Now, in order to measure an atom in the a x basis, we 
simply apply the resonant microwave field to implement 
a rotation of the state by it/ 2 about the y-axis and then 
measure the atom to see whether it is in the state \g) 
or |e). We can see that this works because the rotation 
maps the eigenstates of a x directly onto the states \g) and 
|e). Similarly, to measure in the a y basis we apply a n/2 
rotation about the z-axis (using the Stark shift), then 
a 7r/2 rotation about the y-axis, and then measure the 
atom to see whether it is in \g) or \e) . Measuring Rydberg 
atoms is discussed by Gallagher [21| . The particular mea- 
surement process that is useful in this experiment, state 
selective field ionisation, is described in [17[ . 

So far, we have considered the case of detecting N- 
partite entanglement by making measurements on the 
total system. However, it should also be possible to de- 
tect multiparticle entanglements on subsystems. Suppose 
we consider only the atoms that pass through the first N 
cavities in Figure [T] We ensure that all the a y mea- 
surements are made on these N atoms and that only a x 
measurements are made on the remaining atoms. In this 
case, the expectation value for a general term in (J3j has 
the form, 

First N sites 

= fi{$\4-44 +1 -<£\*)s = i i , (w) 

where any permutation of the first N operators in ei- 
ther expression does not change the result and j is even. 
From ([33]), we see that restricting all the a y measure- 
ments to the first N atoms gives the same results as the 
full Mermin inequality for a state of the form of fl} , but 
with N (rather than N) atoms. From this observation, 
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it follows that we have F = 2 N ~ 
(N — l)-partite entanglement is 



F < 2 



N-2 



and the bound for 



(15) 



The corresponding bound for the hidden variable 
model follows easily from ^ since the joint probabil- 
ity for N sites can be reduced to a joint probability for 
N sites simply by integrating over m^ r+1 , • • • , rnjy. This 
gives, 



P, 



.mi 



'N 



= J d\p(X) [p^mi, A) (77^,^,(16) 

i.e. the first N sites do not depend on the measurement 
outcomes of the last N — N sites. The rest of the ar- 
gument follows that of Mermin [2j but with N replaced 
with N. This gives the hidden variable bounds as, 



F < 2^/ 2 , N even, 

F < 2^" 1 )/ 2 , N odd. 



(17) 



Comparing (fT7|) and (|15|) with the result for a GHZ state, 
F = 2 , we see that it is possible to rule out local- 
realistic models and detect TV-partite entanglements on 
a part of the total state. 

Finally we ought to comment on how noise affects these 
results since it is well-known that detector inefficiencies 
in particular give rise to the so-called detection loophole 
which can undermine our ability to exclude local-realistic 
descriptions using Bell- type inequalities. Braunstein and 
Mann [22[ have considered this problem and shown that if 
the noise is sufficiently small, then the signal for violation 
grows exponentially faster with N than the noise. In 
particular, the noise per detector or per particle needs to 
be less than about 14%. This bodes well for the feasibility 
of the detection process. 

Besides its significance in tests of quantum mechan- 
ics versus local realism, the scheme we have proposed 
could be an important tool in unambiguously creating 
and identifying genuine ./V-body entanglement in atomic 
systems. 
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